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Abstract. Although 70% of the stars in the Galaxy are M-dwarfs, thermal emission searches for cold debris disks have been 
conducted mostly for A-type and solar-type stars. We report on new X = 1.2 mm continuum observations of thirty M-dwarfs, 
using the MAMBO-2 bolometer array camera at the IRAM 30m telescope. For a statistical analysis, we combine these data 
with our prior SCUBA and MAMBO-2 observations of 20 other M-dwarfs. Our total sample divides in M-dwarfs in moving 
groups, with relatively young ages, and in nearby M-dwarfs with unknown ages. Only one cold debris disk (GJ842.2) was 
detected significantly. We compare the implied disk abundance constraints with those found in two comparable submillimeter 
surveys of 10 to 190 Myr old A- and FGK-type stars. For the 19 youngest (ages less than 200 Myr) M-dwarfs in our sample, 
we derive a disk fraction of 5.3^° 5 %, compared to 15+} jj % for FGK-stars and 22 + J 3 ) % for A-stars. Hence, for this age group, 
there is an apparent trend of fewer cold disks for later stellar types, i.e., lower star masses. Although its statistical significance 
is marginal, this trend is strengthened by the deeper sensitivity of observations in the M-dwarf sample. We derive a cold disk 
fraction of < 10 % for the older (likely a few Gyr) M-dwarfs in our sample. Finally, although inconclusively related to a debris 
disk, we present the complex millimeter structure found around the position of the Ml. 5 dwarf GJ526 in our sample. 
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1. Introduction Our current understanding of debris disks has been recently 

reviewed by Wyatt (2008). Mutual collisions between planetes- 

Cold debris disks around main sequence stars are left-over imalfj fa A ^ diskfj produce second . generation dust grains 

£ . planetesimals (comets) that could not agglomerate into larger ^ are observable through their thermal emission or scattered 

planets during the initial phase of planet formation. They are Hght since ^ discovery of a first debris disk the A0 

assembled as a belt in the periphery of the system in a man- main sequence star Vega by IRA S (Aumann et al. 1984), debris 

ner analogous to the Kuiper Belt. The study of debris disks, diskf , haye been searched for photometrically with the infrared 

including warm disks such as the asteroid belt in our Solar (TR) TRAS mQ and Sphzer ^ such obse r Vations , any 

System, advances our knowledge of the origin and evolution flux £xcess abov£ ^ photospheric level is interpreted as emis- 

of planetary systems around other stars, similarly to the study sion from warm (5Q _ m R) circumstel i ar dust . 

of the Kuiper Belt. It is well recognised that its present-day 

structure and dynamics retain important information on the for- Through Spitzer observations it was found that 32 ± 5% 
mation and evolution of the Solar System. For example, the of 160 A-dwarfs show a 70 //m excess (Su et al. 2006), while 
low mass and expanded size of the present-day Kuiper Belt can only 16+|;|% of 225 observed FG-dwarfs show excess emis- 
be traced back to the outward migrations of the giant planets, sion (Bryden et al. 2006; Trilling et al. 2008). Submillimeter 
which exchanged orbital energy with an initially more compact photometry has shown that some A-to-G type stars with no de- 
and more massive disk (Hahn & Malhotra 1999, Tsiganis et al. tectable IRAS excess do however show cold (10 - 50 K) dust 
2005, Morbidelli et al., 2005, and Gomes et al. 2005). emission (Wyatt et al 2003; Najita & Williams 2005). Imaging 
of scattered light with the HST and of thermal continuum emis- 
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2005; Holland et al. 1998; Greaves et al. 2005) ; 22 debris 
disks have been imaged presently Q. Azimuthal structures have 
been detected in a few of these disks and are thought to be 
caused by dust associated with planetesimals trapped in mean 
motion resonance with an orbiting planet (Wyatt 2003, 2006, 
Reche et al. 2008), or by dynamical perturbations from a distant 
stellar companion or passing stars (HD141569A: Augereau & 
Paploizou 2004). 

Although low-mass M-dwarfs are the most populous (70%) 
stars in the Galaxy, they have so far received little attention, 
mostly because their low luminosity makes the thermal emis- 
sion and scattered light from their disks more difficult to detect. 
In a SCUBA survey of young stars of the /3 Pic moving group 
and of the Local Association, Liu et al. (2004) detected the first 
two debris disks around M-dwarfs, AU Mic and GJ182. Using 
SCUBA and MAMBO-2, Lestrade et al. (2006) surveyed 32 
relatively young M-dwarfs of moving groups and newly de- 
tected one disk around the M0. 5 dwarf GJ842.2. Using Spitzer, 
Gautier et al. (2007) surveyed 62 nearby M-dwarfs at 24 //m, 
and subsamples of 41 at 70 fim and of 20 at 160 fim, and found 
no firm detection. AU Mic was also imaged in scattered light, 
revealing an edge-on, structured disk (Liu, 2004 and Krist et 
al., 2005). 

Here we present new MAMBO-2 observations of nearby 
M-dwarfs, and combined them with our previous survey to an- 
alyze a total sample of 50 M-dwarfs (Table 1) in terms of their 
cold debris disks abundance. 

We present the new M-dwarfs surveyed in §2, we describe 
the observations in §3, and results for debris disks but also for 
background sources in the fields of some of the M-dwarfs in §4. 
Finally, we compare the detected fraction of cold debris disks 
around the M-dwarfs of our sample with the fractions for other 
stellar types to discuss whether or not the characteristics of M- 
dwarfs - mass and luminosity lower than solar - impact the 
formation of debris disks around them or their detectability. 

2. Sample of newly observed M-type dwarfs 

To complement our first survey of M-dwarfs, which are in 
moving groups of ages < 600 Myr (Lestrade et al 2006), we 
observed the most nearby M-dwarfs, irrespective of age. We 
selected single M-dwarfs at a distance less than 6 pc and at 
6 > -11°, and added five M-dwarfs binaries and six single 
M-dwarfs between 6 and 10 pc that are in common with the 
Spitzer survey by Gautier et al. (2007). The five binaries in our 
sample are: GJ725 (M3 and M3.5) separated by 15" (53 AU), 
GJ234 (M4.5 and M8), separated by 1" (5 AU), GJ412 (M2 
and M6), separated by 32" ( 1 60 A U), GJ569 (M2 .5 and M8 .5), 
separated by 1" (10 AU), and GJ65 (M5.5 and M5.5), sepa- 
rated by 2" (5.2 AU). These angular separations are so small 
that a single MAMBO-2 map can cover both components. The 
ages of the targeted near-by M-dwarfs are presently unknown. 

3. Observations 

The diameter usually adopted for debris disks is 120 AU, which 
for near-by stars is larger than the IRAM 30-meter telescope 

1 http://astro.berkeley.edu/ kalas/disksite/pages/gallery.html 



beam of 10.7" FWHM at A — 1.2 mm. We have therefore im- 
aged each field with the 117-channel Max-Planck Bolometer 
array (MAMBO-2; Kreysa et al. 1998) of the 30m telescope 
on Pico Veleta, Spain (2900 m). MAMBO-2 has a half-power 
spectral bandwidth from 210 to 290 GHz, with an effective fre- 
quency centered at 250 GHz (1.20 mm) for thermal emission 
spectra. The effective FWHM beam is 10.7", and the under- 
sampled field of view of the array is 4'. We used the standard 
on-the-fly mapping technique, where one map is made of 41 
azimuthal subscans of 60 sec each, with a scanning velocity of 
4"/sec and an elevation incremental step of 4", while chopping 
the secondary mirror at 2 Hz by 60" in azimuth. The bolome- 
ters are arranged in a hexagonal pattern with a beam separation 
of 22". This scanning pattern produces time streams of data 
that are converted to a fully sampled spatial map with 3.5" 
pixels. Our observations were done within pooled observing 
runs spread over the winter and summer periods from 2005 
to 2007. Atmospheric conditions were generally good during 
the observations, with typical zenith opacities between 0. 1 and 
0.3 at 250 GHz and low sky noise. The telescope pointing was 
checked before and after each map by using the same nearest 
bright point source, and was found to be stable to better than 
3", except in a few occasions for which we discarded the data. 
The absolute flux calibration is based on observations of sev- 
eral standard calibration sources, including planets, and on a 
tipping curve (sky dip) measurement of the atmospheric opac- 
ity once every few hours. The resulting absolute flux calibration 
uncertainty is estimated to be about 10% (rms). 

The data were analyzed using the mopsic software pack- 
age written by R. Zylka at IRAM. The double-beam maps 
were combined to a single map using the shift-and-add pro- 
cedure. Compared to a proper image restoration this produces 
maps with about a factor 2 better sensitivity, at the expense 
of no sensitivity to emission structures in scan direction that 
are larger that the wobbler throw of 60". In our sample, the 
shortest integration time per field is 30 minutes, yielding an 
rms noise level of ~ 2.5 mJy/11" beam in the central part of 
the map (r < 60"), steadily rising to ~ 5 mJy/11" beam at 
its edge (r ~ 140"). This non-uniform noise across the maps 
results from the fact that the scanned field is about twice as 
large as the bolometer array size, so that more data is taken 
in the central part of the map than near the edges. Due to this 
non-uniformity, it is judicious to present Signal-to-Noise maps 
rather than intensity maps. The longest observation, by far, had 
20 hour duration and was on GJ628, yielding an rms noise level 
of ~ 0.5 mJy/11" beam in the central part of the map and 
~ 1.2 mJy/11" beam at r ~ 140". The sky area covered by 
each map corresponds to ~ 1700 independent antenna beams. 
The residual noise in the map was found to be nearly Gaussian. 
Hence, a 4<x detection in our maps is statistically significant. 
Below, we describe the procedure we used to extract discrete 
sources and to search for extended emission. 

4. Results 

Although initially intended, our complete survey of 50 M- 
dwarfs turns out not to be flux-limited, which makes its sta- 
tistical interpretation not straightforward. First we present the 
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Table 1. The 50 M-dwarfs of our two (sub)mm surveys. Survey I was already presented by Lestrade et al. (2006) but its data 
were included in the statistical analysis of this paper. New data were acquired in survey II with MAMBO-2. Selection criteria are 
ages < 600 Myr for survey I and the nearest M-dwarfs but of undetermined ages for survey II. Ages are based on Moving Groups 
(Local Association (20 - 150 Myr), IC2391 (35 - 55 Myr), AB Dor (100-125 Myr), Castor (200 Myr), Ursa Maj (500 Myr) and 
Hyades (600 Myr)) identified by Montes et al. (2001) and Zuckerman & Inseok Song (2004a,b). Some stars were observed in 
both surveys. 



Name 


Spectral 


Binarity 


Selection 
criterium 


Observed in 
survey # 


Bolometer 
arrays 


Publication 


GJ82 


M4 


Single 


20- 150 Myr 


I 


SCUBA 


Lestrade et al. (2006) 


GJ212 


M0.5 






I 






GJ507.1 


M1.5 






I 






GJ696 


M0 






I 






GJ876 


M4 






I 


MAMBO 




GJ628 


M3.5 






I 






GJ402 


M4 






I 






GJ234 


M4.5+M8 


Binary 




I & II 




Lestrade et al. (2006) and this work 


GJ285 


M4.5 


Single 




I 


SCUBA & MAMBO 


Lestrade et al. (2006) 


GJ393 


M2 






I 






GJ9809 


M0 






I 






GJ875.1 


M3 




35 - 55 Myr 


I 


SCUBA 




GJ856A 


M3 




100-125 Myr 


I 


MAMBO 




GJ277B 


M3.5 




200 Myr 


I 


SCUBA 




GJ842.2 


M0.5 




?> 


I 




» 


GJ890 


M2 






I 






GJllll 


M6.5 




»» 


I & II 


MAMBO 


Lestrade et al. (2006) and this work 


GJ408 


M2.5 






I & II 






GJ4247 


M4 






I & II 


SCUBA & MAMBO 




GJ447 


M4.0 




500 Myr 


I & II 


MAMBO 




GJ625 


M1.5 






I & II 






GJ569 


M2.5+M8.5 


Binary 




I & II 






GJ873 


M3.5 


Single 




1 




Lestrade et al. (2006) 


GJ65 


M5.5+M5.5 


Binary 


600 Myr 


1 






GJ3379 


M4 


Single 




1 






GJ849 


M3.5 






1 






GJ791.2 


M4.5 






1 






GJ109 


M3.5 






I & II 




Lestrade et al. (2006) and this work 


GJ699 


M4.0 


Single 


1.82 pc 


II 


MAMBO 


this work 


GJ406 


M6.0 




2.38 pc 


II 






GJ411 


M2.0 




2.54 pc 


II 






GJ905 


M5.5 




3.16 pc 


II 






GJ725 


M3+M3.5 


Binary 


3.57 pc 


II 






GJ54.1 


M4.5 


Single 


3.72 pc 


II 






GJ273 


M3.5 




3.79 pc 


II 






GJ83.1 


M4.5 




4.44 pc 


II 






GJ687 


M3.0 




4.53 pc 


II 




» 


LHS292 


M6.5 




4.54 pc 


II 






GJ1002 


M5.5 




4.69 pc 


II 






GJ412 


M2+M6 


Binary 


4.83 pc 


II 






GJ388 


M3.0 


Single 


4.89 pc 


II 






GJ445 


M3.5 




5.38 pc 


II 






LHS1723 


M4.5 




5.43 pc 


II 






GJ526 


M1.5 




5.43 pc 


II 






GJ251 


M3.0 




5.57 pc 


II 






GJ205 


M1.5 




5.71 pc 


II 






GJ213 


M4.0 




5.87 pc 


II 






GJ908 


M1.0 




5.93 pc 


II 






GJ581 


M3 




6.27 pc 


II 






GJ102 


M4 




7.75 pc 


II 
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discrete millimeter sources detected in four MAMBO-2 maps 
(Figs 1, 2 and 3, and Table 2), and we discuss the nature of the 
intriguing cluster of sources around GJ526. Second, we present 
the deep search for faint debris disks made by averaging inten- 
sities over an effective area in the 42 MAMBO-2 maps of our 
complete survey (Tables 3 and 4). The 8 other M-dwarfs were 
observed in wide photometry with SCUBA (Table 5) and were 
already discussed by Lestrade et al. (2006). We use the com- 
plete sample of 50 M-dwarfs to estimate the fraction of cold de- 
bris disks around M-dwarfs and upper limits of their fractional 
dust luminosities. Three stars (GJ285, GJ393 and GJ4247) are 
in common between Tables 3, 4 and 5. 

4.1. Discrete sources and structure around GJ526. 
4.1 .1 . Source extraction 

In Table 2, we summarize the characteristics of the discrete 
sources detected with S /N > 4 in the 42 fields mapped with 
MAMBO-2, each being ~ 400" x 400" in size and centered on 
the position of an M-dwarf. The source extraction was done by 
searching each map for any pixel with S /N > 4 to select a block 
of 7 X 7 pixels centered on it. Then, we minimized^- 2 between 
the 2-D Gaussian F x exp( - 0.5 x [(* - x ) 2 + (y- yo) 2 ]/o" 2 ) 
and the measured intensities over this block by varying the peak 
flux F, the parameter cr (FWHM= 2.35<x) and the peak posi- 
tion (xo,yo) by less than a pixel from the S /N > 4 pixel. We 
used the mean brightness computed from the 32 pixels sur- 
rounding the 7x7 block to remove any constant level. With 
this scheme, this level is determined from pixels that are far 
enough from the 2-D Gaussian peak to be at the map floor 
since the Gaussian FWHM is 3 to 5 pixels (11" - 17") for 
our sources. This correction amounted to between 5% and 20% 
of the integrated flux density. The integrated flux densities of 
sources in Table 2 are S „, = 2ncrF '. We have extracted 13 dis- 
crete sources with S /N > 4 for F in 4 fields out of the 42 
MAMBO-2 maps ; 8 sources are newly found in the fields of 
GJ526, GJ725 and GJ569. They are shown in Fig. 1, 2, and 
3. The 5 others sources are in the field of GJ628 as already 
reported in Lestrade et al. (2006) ; typographic errors in their 
coordinates in this first publication are corrected in Table 2. In 
this Table, the lowest integrated flux density is 3.0±0.7 mJy and 
the highest is between 30 and 63 mJy. There is no optical coun- 
terpart within 7" (2<x MAMBO-2 position error) to any of these 
sources in the USNO-B1 catalogue (Monet et al. 2003), except 
for MM 163007- 123942 (GJ628W) at the lcr level in position 
(object I mag = 18.92) as already reported in Lestrade et al. 
(2006). Additionally, we have searched the fields of GJ526 and 
GJ569 that have been imaged in the Sloan Digital Sky Survey 
(SDSS). In the field of GJ526, there is only one SDSS object (z 
mag = 20.4) which is as close as 7" from MM134540+145446. 
However, the confusion limit is reached over such a separa- 
tion since the density of SDSS objets is -22000 objects/deg 2 
as measured over the 4' x 4' area centered on the star position. 
In the field of GJ569 (-33000 SDSS objects/deg 2 ), there are 
two SDSS objets (mag u =19.9, g=21.6) as close as 8" from 
MM 145428 +160439, but again confusion prevents identifica- 
tion with the MAMBO-2 millimeter source. We have found 



no radio counterpart within the 3cr MAMBO-2 position un- 
certainty (10.5") in the NVSS catalogue (Condon et al., 1998) 
to the flux density limit of 2.5 mJy at 1.4 GHz, except for 
MM163007- 123942 (GJ628W) at the 1.7<x level in position 
(object flux density = 5.6 ± 0.5 mJy) as already reported in 
Lestrade et al. (2006). Finally, we have found no IR counterpart 
within the 3<x MAMBO-2 position uncertainty in the 2MASS 
All-Sky Catalog of Point Sources (Skrutskie et al., 2006). 

4.1 .2. A bright millimeter source 

The source MM 184222+593828 in the field of GJ725 has a flux 
density between 30 and 63 mJy (Table 2), which is bright com- 
pared with millimeter sources found previously in empty fields 
mapped with MAMBO-2 (Greve et al., 2004, Bertoldi et al., 
2007 and Voss et al. 2006). These lower and upper limits were 
obtained by using different atmospheric correlation lengths in 
the skynoise reduction within the mopsic data reduction pack- 
age, yielding similar detections ofS/N~ 7. Systematics affects 
the accuracy of the flux measurement for this source, which is 
near the edge of the map. This source is possibly resolved with 
a source FWHM between 12" and 15" depending on the re- 
duction parameter used. Although this is rare, it has recently 
been shown that submm galaxies can have multiple compo- 
nents (Tacconi et al., 2008). We shall discuss this source to- 
gether with complementary observations in a forthcoming pa- 
per. 

4.1.3. Structure around GJ526 

The map of GJ526 is shown in Fig. 1. Five sources are de- 
tected with S IN > 4 and an ellipse almost centered on the 
star can connect all but the one to the SW. Precisely, there are 
four robust sources with S IN > 4 that can be connected, and 
an additional one, farthest to the SE that is only detected at 
3.5 times the rms noise level, and possibly real since statis- 
tically less than 1 positive noise peak is expected at this level 
over the whole map if dominated by gaussian noise. To estimate 
whether this is the case or not, we have identified the negative 
noise spikes S /N < -3 and the positive spikes 3 < S /N < 4 
in our 440" x 420" map (~ 2050 beams). The total of 7 such 
noise spikes found in the map is consistent with the 6 theo- 
retically expected for \S /N\ > 3 for gaussian noise. Also, the 
presence of the SE source is found robust since it appears at 
the same position in the two independent maps we made by 
halving the data. There is also a ~ 3cr source to the West close 
to the ellipse but this source is not robust to data selection. So 
we have disregarded it. 

The ellipse is offset by two pixels in right-ascension from 
the position of the star, its semi-major and semi-minor axes 
are a = 98" and b = 15", and its orientation is PA=-30°. 
This structure might be the projected ring of a clumpy debris 
disk, its inclination being 81° from the plane of the sky. The 
offset between the star and the center of the ellipse is only 2 
pixels and might be caused by position uncertainties or/and by 
real source structure. The 5 connected sources appear not to 
be embedded in any extended emission as expected for a de- 
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Table 2. Sources found in our MAMBO-2 survey (S/N > 4). The IRAM30m beam is 11" at A = 1.2 mm. The flux density of 
MM184222+593828 is provided as lower and upper limits because the source is located close to one of the borders of the map 
where some systematics remain in the data (see text). 



OULllLC lltllllL 


Star Field 

..J l ill 1 1C1U 


a(J2000) " 




TntpOTiitprl flii"Y HptiQitv ^ 

llllCgl cllCU 1 1 LI A UG11B11Y 

nf 1 9 mm (mhn 
at l.Z, 111111 ^llljy^ 


S IN 


Qniirpp 

. ) KJ HI V- C 

r VV JTLIVI 


MM145428+160439 


GJ569 


14 


54 


28.2 


16 04 39 


5.3 ±1.1 


4.9 


ii" 


MMI SzL999-i-SQ^S9S 
1V11V1 1 OtZZZt J^JJOZo 


^ J.I / LJ 


is 


42 


IL. 1 


38 98 
Jy jo ZO 




fx 1 7 3 


1 V 1 V 
lj — 1 j 


IVllVl 1 OtZJJtjyj / JO 


OT79S 


18 


42 


53.3 


SO 37 


Q 7 + 1 Q 


S 1 

J. 1 


1 1 


MM134539+145139 


GJ526 


13 


45 


39.0 


14 51 39 


3.2 ±0.7 


4.4 


ii" 


MM134540+145446 


GJ526 


13 


45 


40.7 


14 54 46 


4.3 ± 0.8 


5.2 


ii" 


MM134541 + 145417 


GJ526 


13 


45 


41.3 


14 54 17 


6.3 ±1.0 


6.1 


ii" 


MM134543+145317 


GJ526 


13 


45 


43.1 


14 53 17 


3.0 ±0.7 


4.9 


ii" 


MM134546+145240 


GJ526 


13 


45 


46.7 


14 52 40 


5.6 ±0.7 


8.6 


ii" 


MM163007-123942 


GJ628 c 


16 


30 


07.6 


-12 39 42 


7.3 ± 1.2 


6.2 


14 ± 4" 


MM163022-123925 


GJ628 d 


16 


30 


22.3 


-12 39 25 


3.8 ±0.9 


4.4 


17 ± 5" 


MM163019-123830 


GJ628 e 


16 


30 


19.7 


-12 38 30 


4.1 ±0.8 


4.9 


11" 


MM163015-123911 


GJ628 f 


16 


30 


15.6 


-12 39 11 


4.8 ±0.8 


6.2 


11" 


MM163013-124057 


GJ628 s 


16 


30 


13.7 


-12 40 57 


4.3 ± 1.0 


4.8 


11" 



" The MAMBO coordinate uncertainties are ~ 3.5". 

h Flux density (mJy) is integrated under the fitted 2D-Gaussian. 

c , d , f , 8 : these sources were found in survey I ; ( c ) is GJ628-W, ( d ) is GJ628-E, (') is GJ628-NE, (f) is GJ628-NW, and (*) is GJ628-SW in Fig 2 and Table 3 of Lestrade et al. 
(2006). Note that there are typographic errors in the minute column of declinations in Table 3 of Lestrade et al. (2006). The positions of the GJ628 sources given now are correct and 
supersede this first publication. Also, the flux densities in this first publication were peak flux densities while they are integrated flux densities now. 



bris disk though. The mean brightness over a rectangular box 
220" x 40" (92 beams) oriented at PA =-30° and centered on 
the star is 0.21 ± 0.073 mJy/1 1" beam, i.e. ~ 3cr. Actually, this 
mean brightness matches the mean of the 5 flux densities of the 
connected sources within the box, indicating that any extended 
emission must have a brightness < 3 x 0.073mJy/ll" beam at 
1.2 mm. 

We now test whether or not such a source cluster can arise 
from the distribution of background sources. We carry out a 
first test to estimate the probability that the connected sources 
around GJ526 can cluster as tightly as they do in Fig 1 if they 
were background sources. For this test, we use the statistical 
analysis of spatial point patterns developed by Diggle (2003). 
This analysis is based on the nearest neighbour distance, de- 
fined as the distance between a point (a source for us) and its 
nearest neighbour. For GJ526, all the nearest neighbour dis- 
tances for the 5 sources connected by the ellipse are < 60". The 
probability that k distances x be < 60" among distances if the 
mean frequence of occurence for x < 60" over the whole sky is 

/, is given by the binomial distribution B k = ( ) J* (1 - f) N ~ k . 
We can derive the mean frequence of occurence / from the 
three empty fields (ELAIS N2, Lockman Hole and COSMOS) 
mapped by MAMBO-2 (Greves et al., 2004 and Bertoldi et al., 
2007). We found there are 9 distances < 60" between the 71 
sources of these three fields, and so / = 9/(71 - 3) = 0.13. 
The factor (71-3) is the number of nearest neighbour distances 
for 71 sources distributed in the 3 separate fields. For the five 
sources around GJ526 in Fig 1, all 5-1 nearest neighbour dis- 



tances are < 60" as already mentioned, so k = N = 4 and 
the probability that background sources can produce such a 
number of small distances is as low as 0.1%. 

As a second test, we use the elliptical pattern connecting 
6 submm sources recognisable in the North- West part of the 
COSMOS field mapped by MAMBO-2 and displayed in Fig 2 
of Bertoldi et al. (2007). This field is 20x20 arcmin 2 in size 
and the mean occurence of ellipse of major-axis 200", as in the 
field of GJ526, is 1 /20x20 arcmin -2 . We compute the Poisson 
probability to find a similar ellipse in a field that we take as 
small as 200x200 arcsec 2 to account for the fact that the ellipse 
around GJ526 is found in an area restricted to the central part of 
the map, i.e. centered on the position of GJ525. This probability 
is 2.7%. This test is only indicative because the COSMOS field 
has been mapped at the level of 1 mJy/beam rms while our map 
is twice as deep for the field around GJ526. 

In summary, the five sources in Fig 1 that are symmetricly 
located around GJ526 are connected by an ellipse almost cen- 
tered on the star. The two tests carried out above provide indica- 
tions that this structure is statistically unconsistent with known 
spatial distributions of background submm galaxies. Hence, at 
this stage, we cannot rule out the hypothesis that the sources are 
associated with the star. In this case, the 5 connected sources 
could be indicative of azimuthal structures in an inclined de- 
bris disk around GJ526 whose extended emission is not seen 
because the map is not deep enough. Complementary obser- 
vations at 850 //m, and at shorter wavelengths with Herschel 
should attempt to detect the extended emission of the disk. 
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Table 3. New MAMBO-2 observations at 1.2 mm for nearby M-dwarfs. The table includes mean brightness determined by 
averaging map intensities over an effective area to search for faint debris disks that are not readily apparent. Not shown in this 
table are the five emission clumps symmetricly located around GJ526 that might be an inclined debris disk (see subsection 4.1.3) 
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M1.5 


5.71 


-19 


1.0 
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< 2.0 


GJ213 


M4.0 


5.87 


-9 


0.5 


2.19 


-0.31 ±0.44 


-0.70 


< 2.4 
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M1.0 


5.93 


-57 


0.5 


2.01 


0.24 ± 0.39 


0.61 


<2.2 


GJ581 


M3 


6.27 


40 


0.5 


1.96 


-0.75 ± 0.43 


-1.75 


<2.1 


GJ625 


M1.5 


6.58 


43 


5.0 


0.69 


0.38 ±0.15 


2.53 


<0.7 


GJ408 


M2.5 


6.62 


64 


2.0 


1.13 


0.42 ±0.14 


3.05 


< 1.2 


GJ109 


M3.5 


7.55 


-31 


3.0 


0.98 


0.32 ± 0.22 


1.48 


< 1.1 


GJ102 


M4 


7.75 


-32 


0.5 


2.20 


-0.09 ± 0.60 


-0.15 


< 2.4 


GJ4247 


M4 


8.96 


-21 


3.5 


0.74 


0.63 ± 0.25 


2.56 


<0.8 


GJ569 (*) 


M2.5+M8.5 


9.81 


59 


2.0 


1.09 


-0.46 ± 0.37 


-1.25 


< 1.2 



(1) rms estimated for r < 60" in the map ; 

(2) mean brightness and uncertainty of mean computed by averaging intensities over an effective disk of radius 60 AU ; 

(3) the 3o" flux density upper limit is computed also over the same effective area (radius=60AU). 
(*) : background sources in map (see Table 2) ; 

(**) : possibly a large debris disk (see Fig 1 and subsection 4. 1.3). 



Eventually, astrometry should detect the same prope motion 
for the 5 connected sources as for the star GJ526 (2.3 "/yr in 
the SE direction) if indeed they are part of a debris disk. The 
mid-epoch of our MAMBO-2 data is early 2007, so that the 
5 sources should have moved in concert with the star by a full 
IRAM-30m beam by 20 1 1 , providing definitive proof of a disk. 

Such a debris disk would have a radius as large as ~ 500 AU 
at the distance of GJ526. We examine whether or not this 
is conceivable. First, we note that exceptionally large debris 
disks, 520 AU and 600 AU in extent, have recently been found 
around the 184 Myr old AO-type dwarf y Oph at 70 //m (Su et 
al. 2008) and around the main sequence F8V q l Eri at 870 fim 
(Liseau et al. 2008). Second, protoplanetary disks where plan- 



etesimals and planets form extend to almost 1000 AU as ob- 
served for example around the young close binary GG Tau 
(Dutrey, Guilloteau & Simon, 1994). In the model proposed by 
Kenyon and Bromley (2004a), icy-planets successively form in 
waves outward in the disk producing larger and larger dusty 
rings from collisional cascades. In their model, the planet for- 
mation timescale is 15 - 20 x (2 /2 MMSA ,) _1 (a/30 AU) 3 Myr 
in a quiet disk (their eq (4)). Assuming for Eo as much as 15 
times the surface density of the minimum-mass solar nebula 
Z-mms n as required for the formation of Jupiter in the solar sys- 
tem (Lissauer 1987), the timescale for planet formation to reach 
500 AU in the GJ526 system is ~ 5 Gyr. The stellar diame- 
ter of GJ526 measured by the optical interferometer CHARA 
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Table 4. M-dwarfs associated with moving groups observed at 1.2 mm in survey I with the MAMBO-2 facility and already 
presented in Lestrade et al. (2006) but reanalyzed here in a consistent fashion with the new data of Table 3. 



Name 


Spectral 


D 


Gal. 1. 


Integration 


Map rms 


Mean brightness 


u/o~,, 


3 (X flux density 




tvDe 


(vc) 




time 


(mJy/11" beam) 


LL + o~,, 


(2) 


upper limit (1.2 mm) 










(hours) 




(mJy/11" beam) 




(mJv) 












(1) 


(2) 




(3) 


GJ65 


M5.5 


2.6 


-76 


1.0 


2.30 


-0.15 ± 0.17 


-0.87 


< 10.4 


GJ628 


M3.5 


4.5 


24 


20.0 


0.51 


-0.06 ± 0.08 


-0.74 


< 1.3 


GJ876 


M4 


4.7 


-59 


0.5 


2.80 


0.17 ±0.46 


0.37 


<7.0 


GJ873 


M3.5 


5.0 


-13 


0.5 


2.31 


0.33 + 0.38 


0.86 


<5.4 


GJ3379 


M4 


5.4 


-10 


0.5 


2.57 


0.05 ± 0.56 


0.09 


< 5.6 


GJ285 


M4.5 


5.9 


13 


1.0 


1.74 


-0.88 ± 0.35 


-2.52 


< 1.9 


GJ402 


M4 


6.8 


55 


1.0 


1.81 


-0.86 ± 0.39 


-2.18 


< 2.0 


GJ393 


M2 


7.2 


47 


1.5 


1.51 


-0.50 ± 0.33 


-1.50 


< 1.6 


GJ849 


M3.5 


8.8 


-45 


0.5 


2.71 


0.66 ± 0.89 


0.74 


<2.9 


GJ791.2 


M4.5 


8.9 


-17 


1.0 


1.70 


-0.87 ± 0.57 


-1.53 


< 1.8 


GJ856A 


M3 


16.0 


-20 


1.1 


1.50 


-0.56 ± 1.50 


-0.37 


< 1.6 


GJ9809 


MO 


24.0 


03 


0.5 


2.20 


0.36 ± 2.20 


0.16 


<2.4 



(1) rms estimated for r < 60" in the map ; 

(2) mean brightness and uncertainty of mean computed by averaging intensities over an effective disk of radius 60 AU ; 

(3) the 3<x flux density upper limit is computed also over the same effective area (radius=60AU). 



(Berger et al., 2006) is about 2<x larger than the ZAMS diam- 
eter predicted by Chabrier & Baraffe (1997) and Siess et al. 
(2000). If such a deviation of only 2cr is real, it is an indica- 
tion of youth instead for GJ526, and this would pose a problem 
for the Kenyon and Bromley model applied to the hypothetical 
debris disk around GJ526. Note that, with the dust surface den- 
sity profile 15I.mmsn x (r/ro) 3 ^ 2 , the mass of solids in a ring 
at r = 500 AU and O.lr in width is 1 M ffi which is enough to 
form planetesimals. Finally, the dust mass corresponding to the 
emission of the 5 sources is between ~ 6 and ~ 10 lunar masses 
as derived in the Appendix for dust grains with a Dohnanyi size 
distribution and heated both by the stellar luminosity and the 
interstellar radiation field. The dust mass is ~ 22 lunar masses 
if derived conventionally for grey body dust with a mass opac- 
ity of 1.7cm 2 g _1 at 850 fim and oc A' 1 for A > 210 fim, and 
with the single dust temperature 4.9K at large radius from the 
star (r = 500 AU) where the interstellar radiation field domi- 
nates the heating process of grains (see Appendix). 



4.2. Deep search for debris disks in the MAMBO-2 
maps 

We searched for faint debris disks in each MAMBO-2 map by 
averaging intensities over a disk of increasing radius till 30" 
to see whether or not the mean brigthness peaks at some an- 
gular radius 6. The radius limit of 30" comes from the shift- 
and-add reduction method and the wobbler throw used of 60" 
for the observations. The intention with this averaging was to 
find a disk whose structure is not directly apparent in the map 
but whose mean brightness is statistically significant. Naturally, 
no information on its structure can be recovered with this pro- 



Table 5. M-dwarfs observed at 850 //m with JCMT/SCUBA 
and already presented in Lestrade et al. (2006) but reanalyzed 
here in a consistent fashion with the new data of Tables 3 and 
4. 



Star 


Sp. 


Dist. 


Integration 


Flux density " 


Size* 




type 


(pc) 


time 


850 yum 


(AU) 








(hrs) 


(mJy) 




GJ82 


M4 


12.0 


2 


2.0 ± 1.4 


84 


GI212 


M0.5 


12.5 


2 


1.3 ± 1.4 


88 


GI285 C 


M4.5 


5.9 


1 


-0.7 ± 1.9 


41 


GJ393 C 


M2 


7.2 


1 


1.9+ 1.9 


50 


GJ507.1 


M1.5 


17.4 


1 


-0.4 ± 2.0 


121 


GI696 


M0 


21.9 


2 


0.8+0.8 


150 


GJ9809 C 


M0 


24.9 


1 


-5.2 ± 2.3 


174 


GJ4247 C 


M4 


9.0 


1 


1.1 ± 2.1 


63 


GJ277B 


M3.5 


11.5 


1 


-2.1 ± 1.8 


80 


GJ842.2 


M0.5 


20.9 


4 


25.0 + 4.6 


300 


GI890 


M2 


21.9 


2 


-2.6+ 1.6 


153 


GJ875.1 


M3 


14.2 


1 


0.52 + 2.1 


99 



1 Flux density from wide photometry. 

' Size is the radius of a debris disk as large as the telescope beam (14"). 
Stars also observed by MAMBO-2 in Tables 3 and 4. 



cedure. The optimum sensitivity of this method is for face- 
on disks while highly inclined disks might escape detections. 
This method is similar to the one used to determine extension 
limits of debris disks in the mid-IR surveys of Sun-like stars 
conducted by Smith, Wyatt & Dent (2008), although, we do 
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Fig. 1. MAMBO-2 Signal-to-Noise ratio map of the field 
around the Ml. 5 dwarf GJ526 at A = 1.2 mm. The pixel size 
is 3.5" x 3.5". The noise rms is <x = 0.58 mJy/11" beam 
in the central region (r < 60") and increases towards the 
edges of the map, ~ 1.3 mJy/11" beam at r ~ 140". 
The contours are —Act, — 3cr, -2<x, -lcr (dotted lines), and 
lcr,2cr,3cr,4cr,5cr,6cr,7cr,8cr. The ellipse (a = 98", b = 15" 
and PA=-30°) is almost centered on the star position and con- 
nects five sources that might be clumps of a debris disk inclined 
to the plane of the sky by 81° (see subsection 4. 1.3). The field is 
centered on the position of GJ526 of early 2006 : a(J2000) = 
13h 45m 44.52s and 5(72000) = 14°53'20.6" (red star). See 
Table 2 for integrated flux densities and coordinates of sources. 



not have the complication of having to accurately substract the 
photosphere at A — 1.2 mm. Following this procedure, we plot- 
ted the mean brightness as a function of angular radius 9 for 
each star of Tables 3 and 4. All the curves were inspected and 
found to wander around zero mean with excursions < 3<x for 9 
comprised between 1 1"/2 and 30", indicating that no new disk 
was found by this method. 

In Fig 4, we provide the distribution of the mean bright- 
ness /j I uncertainty of mean cr^ listed in Tables 3 and 4 and 
computed for 9 corresponding specificly to the adopted disk 
radius 60 AU at the distances of the stars. The corresponding 
Gaussian probability density function is also plotted in Fig. 4. 
Comparison between the two distributions indicates there are 
more high positive and negative ratios /v/cry, than expected al- 
though they stay within ±3<x. It means that there are still some 
systematic errors in the maps but at a low level, likely caused 
by remaining atmospheric fluctuations. 




200 100 -100 -200 

Aaxcos(6) offset (") 



Fig. 2. MAMBO-2 Signal-to-Noise ratio map of the field 
around the binary GJ725 (M3 + M3.5 dwarfs) at A = 
1.2 mm. The pixel size is 3.5" x 3.5". The noise rms is cr — 
1 .4 mJy /ll" beam in the central region (r < 60") and increases 
towards the edges of the map, ~ 4.3 mJy/11" beam at r ~ 
170". The contours are — Act, -3<t, -lcr, -lcr (dotted lines), and 

lcr,2<T,3o-,4o-,5o-,6o-,7<x. The source MM184222+593828 
westward has an integrated flux density between 30 and 63 mJy 
(S/N ~ 7). The other source, MM184253+593756,has an inte- 
grated flux density of 9.7mJy and a S /N of 5.1. The field is cen- 
tered on the position of GJ725 A of early 2006 : a(72000) = 18h 
42m 45.48s and 5(72000) =59°38'01. 7" (red star). See Table 2 
for integrated flux densities and coordinates of the sources. 

4.3. Upper limits on fractional dust luminosities and on 
dust masses of the debris disks 

The fractional dust luminosity is the fraction of the stellar radi- 
ation absorbed and reprocessed to the infrared and (sub)mm by 
the dust grains ; it is proportional to the fraction of the sky cov- 
ered by dust as seen from the star (Dominik and Decin 2003). 
We used the Stefan-Boltzmann law - black body emission - to 
estimate dust luminosity, and modified it by emissivity 1/X A , 
with X A = 1 for A < 210 yum and X A = A/210 for A > 210 /jm 
(Wyatt 2008). In these conditions, fractional dust luminosity 
is : 

W = 9xlOV s, kT , SeJT} 
U 2nhv 3 RlTt 

the normalization coefficient is such that the measured flux den- 
sity S g is in Jy, the star distance d in pc, the stellar radius R„ 
in m, the dust and stellar effective temperatures Tj and 7» in K, 
the Planck and Boltzmann constants /; and k in 7 x s and J/K, 
the speed of light c in m/s, and the frequency of observation v 
in Hz. 

The standard argument used to fix dust temperature in 
mid-IR surveys of debris disks is that observations are most 
sensitive to the dust emission peak, and Tj is derived from the 
Wien law. In the (sub)mm range, this law is not appropriate for 
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200 100 -100 -200 



Aaxcos(6) offset (") 

Fig. 3. MAMBO-2 Signal-to-Noise ratio map of the field 
around the binary GJ569 (M2.5 + M8.5 dwarfs) at A = 
1.2 mm. The pixel size is 3.5" x 3.5". The noise rms is cr = 
1.1 mJy/11" beam in the central region (r < 60") and in- 
creases towards the edges of the map, ~ 2.8 mJy/11" beam 
at r ~ 140". The contours are -Act, -3cr, — 2cr, -lcr (dotted 
lines), and lcr, 2cr, 3cr, 4<x. The source MM 145428 +160439 has 
an integrated flux density of 5.3 mJy and a S /N of 4.9. The 
field is centered on the position of GJ569A of early 2006 : 
or(72000) = 14h 54m 29.35s and 5(72000) = 16°06'03.2" (red 
star). See Table 2 for integrated flux density and coordinates of 
the source. 

Td since we observe in the Rayleigh-Jeans limit. To keep full 
generality, we avoid choosing a dust temperature T c / at some ar- 
bitrary radius but plot in Figs 5 and 6 fractional dust luminosi- 
ties and dust masses as functions of disk radius r comprised 
between 1 and 1000 AU following the approach by Bryden et 
al. (2006) and Wyatt (2008). At disk radius r, we use the dust 
temperature : 

T d = 278 x (L° 25 ) x (r-°' 5 ) (2) 

from black body equilibrium where L* is the stellar luminos- 
ity in L Q and r is in AU (Backman & Paresce 1993). Note that 
by combining eqs (1) and (2), one gets the expression for frac- 
tional dust luminosity given by eq (8) of Wyatt (2008) where 
the stellar luminosity L» cancels out. 

We adopt the 3cr flux density limit for Sg in eq (1) by in- 
tegrating brightness over a face-on disk of radius r, or at the 
distance of the star. If 20 > 11" (i.e.> IRAM 30m beam) : 

Sg = 3 x rms x (20"/ 1 1") (mJy) 

This formula takes into account that the mean brightness un- 
certainty cry, improves as rms/ -^number of beams while the 
integrated flux density increases as number of beams of the 
disk area. Now if 20 < 11" : 

Sg = 3 x rms (mJy) 
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Fig. 4. Distribution of the ratios mean brightness fi / uncertainty 
of mean cr^ for the 42 M-dwarfs observed with MAMBO-2, 
i.e. fi/cr^ listed in Tables 3 and 4. Here, the mean brightness 
fi is computed over a disk of radius 60 AU adopted for each 
star. The Gaussian probability density function superimposed 
is scaled so that its integral is 42. Comparison between the two 
distributions indicates there are more high positive and nega- 
tive ratios yu/o), than expected although they stay within ±3cr. 
It means that there are still some systematic errors in the maps 
but at a low level, likely caused by remaining atmospheric fluc- 
tuations. No mean brightness is retained as statistically signifi- 
cant in the sample. 

We have used the rms of Tables 3 and 4 that correspond to the 
central part of the maps, thus underestimating slightly the upper 
limits computed. 

In Fig. 5, we show the resulting fractional dust luminosity 
upper limits for r comprised between 1 and 1000 AU. These 
functions first show a steep negative slope as long as 20 < 11" 
making Sg constant in eq (1), then these functions level off 
when Sg linearly increases with 0, finally they increase when 
the dust temperature saturates at 4.9K because the interstel- 
lar radiation field becomes dominant over the stellar field (see 
Appendix). In this figure, we have added the upper limits of 
the fractional dust luminosities for the 41 M-dwarfs observed 
at 70 jjm by Spitzer (Gautier et al. 2007), computed in a simi- 
lar fashion from their 3cr flux densities. The figure shows that 
the two sets of data are complementary, and lead to a uniform 
fractional dust luminosity over a large extent of disk radii for 
the dozen of M-dwarfs common to the two data sets. In Fig. 6, 
we present the corresponding dust masses as a function of r 
for both samples computed with the optically thin emission 
model {e.g. Zuckerman 2001) and the mass opacity kh oc A 
for A > 210 /vm and #850/im = 1.7cm 2 g _1 . 

5. Discussion 

The single cold debris disk found in our surveys (GJ842.2, see 
details in Lestrade et al. (2006), and excluding GJ526 at this 
stage) makes the detection fraction to amount to 2+fj % in our 
sample of 50 M-dwarfs. The limits are based on the Binomial 
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Fig. 5. Constraints for the dust luminosity fraction versus disk 
radius. Note that the x-axis shows the single radius correspond- 
ing to the single temperature of a disk (ring) assumed infinitely 
narrow in our model. We consider the wide range of radii from 
1 to 1000 AU as plausible for rings of debris. Upper limits are 
shown as dashed lines, detections as solid lines. Blue lines cor- 
respond to the (sub)millimeter observations ; i.e. our sample of 
50 M-dwarfs for which dark blue is used for "young" M-dwarfs 
(ages < 200 Myr) and light blue for "old" M-dwarfs (likely a 
few Gyr) ; the two submm disk detected around AU Mic and 
resolved (Liu et al., 2004, Liu 2004) marked by a single blue 
dot ; the submm disks detected but not clearly resolved around 
GJ182 (Liu et al. 2004) and around GJ842.2 (Lestrade et al. 
2006) ; the submm transition disk detected around the pre-main 
sequence Ml dwarf TWA7 (Matthews, Kalas, & Wyatt, 2007) 
was also included. Orange lines correspond to Spitzer 70 yum 
observations ; i.e. the sample of 41 M-dwarfs of Gautier et al. 
(2007) and the detection of the M0 dwarf HD95650 (Smith et 
al. 2006). There are 16 M-dwarfs in common between the two 
samples. These two sets of data are complementary to constrain 
the existence of warm dust around M-dwarfs at moderate radii 
(< 20 AU) and cold dust at large radii (> 20 AU). Disk diame- 
ters probed by the observations are limited by the angular size 
of 60" for the MAMBO-2 maps, of 28" for SCUBA wide pho- 
tometry and of 38.4" for Spitzer aperture photometry (Gautier 
et al. 2007). Some curves are terminated at less than 1000 AU 
because of these angular limits. The calculation of the frac- 
tional dust luminosity is described in the text. (This figure is 
available in color in electronic form). 



distribution for a small number sample and are such that 68% 
of the probability is between the lower and upper uncertainties 
and the peak probablity is the observed fraction 1/50, following 
Burgasser et al. (2003). We recall that Gautier et al. (2007) had 
no detection in a sample of 62 nearby M-dwarfs at 24jum and no 
detection in a subsample of 41 of them at 70 yum, i.e. rate < 7%. 
Differently, Forbrich et al. (2008) have detected photometric 
excesses at 24//m tracing warm dust around 1 1 M-dwarfs in 
the young open cluster NGC2457 (20-40 Myr) that represent 
4.9^jg % of the 225 highly probable member M-stars identified 
in it. 



- 1 10 100 1000 

Disk radius (AU) 

Fig. 6. Constraints for the dust mass versus disk radius. Note 
that the x-axis shows the single radius corresponding to the 
single temperature of a disk (ring) assumed infinitely narrow in 
our model. We consider the wide range of radii from 1 to 1000 
AU as plausible for rings of debris. Upper limits are shown as 
dashed lines, detections as solid lines. Blue lines correspond to 
the (sub)millimeter observations and orange lines to the Spitzer 
70 /urn observations as detailed in the legend of Fig. 5. The mass 
opacity used to convert flux density to dust mass with the stan- 
dard optically thin emission model is k = 1.7cmV' at 850 jum 
and oc /T 1 for A > 210 jt/m. Additional information are in the 
text. (This figure is available in color in electronic form). 
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Fig. 7. Detection rates of cold debris disks versus stellar 
masses (stellar types) for stars younger than 200 Myr. These 
rates and uncertainties are for disks having dust fractional lu- 
minosities larger than the limits shown in Fig. 8 of the (sub)mm 
surveys used. 

To discuss how detectability of cold debris disks depends 
on the mass of the central star, we have compared our re- 
sult with the observed fractions of cold debris disks around 
stars more massive than M-dwarfs. In the literature, we find 
~ 30 cold debris disks around A-to-K type stars detected by 
submm observations (see Fig 3 and caption of Wyatt 2008). 
But most of these detections come from JCMT/SCUBA sur- 
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Fig. 8. Dust luminosity versus disk radius limits of the three 
(sub)mm surveys used to determine the cold disk fractions of 
"young" M-, FGK- and A-type stars. Shown are the dust lu- 
minosity fractions of the 19 youngest M-dwarfs of our sam- 
ple (blue), of the 9 young A-stars (green) and of the 26 young 
FGK-stars (red) in the two 850/im surveys of Wyatt, Dent & 
Greaves (2003) and of Najita & Williams (2005). All these 
stars are less than 200 Myr old. Most of the curves are upper 
limits (dashed lines). The 7 detections are marked by full lines. 
This plot emphasizes that for stars of lower masses (later stellar 
types), even though the surveys are more sensitive, fewer disks 
are detected. Note that the x-axis shows the single radius cor- 
responding to the single temperature of a disk (ring) assumed 
infinitely narrow in our model. We consider the wide range of 
radii from 1 to 1000 AU as plausible for rings of debris. (This 
figure is available in color in electronic form). 

veys of IRAS biased samples, i.e. targets with prior IRAS ex- 
cess detections, unlike our M-dwarfs, for which no such prior 
knowledge was used in the selection of the sample. There are 
only two submm surveys of A-to-K type stars that are unbiased 
in this respect and that have depths of ~l-3 mJy at 850 //m, 
comparable to our sensitivity at 1 .2 mm : the JCMT survey by 
Najita & Williams (2005) of thirteen F5-to-K3 stars (10 Myr 
< ages < 180 Myr, 10 < d < 78 pc, 3 detections) and the 
JCMT survey by Wyatt, Dent & Greaves (2003) of nine B7- 
to-AO stars (86 < d < 938 pc, 2 detections) and of thirteen 
F3-to-K5 stars (27 < d < 250 pc, 1 detection) that are all part 
of Lindroos binaries (14 Myr < ages < 170 Myr). Combining 
these two submm surveys, we determine the cold disk fractions 
of 22+g % for young A stars (2 detections/9) and of 15+ J }| % 
for young FGK stars (4/26). For a meaningful comparison, in- 
structed by Spitzer surveys where the age factor is crucial for 
detectability of warm debris disks, we selected the 19 youngest 
M dwarfs which have ages between 20 and 200 Myr in our sam- 
ple. The disk fraction for these "young" M-dwarfs is 5.3^° 5 %. 
Therefore, for this age range, there is an apparent trend in these 
three fractions, indicative of fewer cold disks detected for later 
stellar types - lower star masses - although at a low statistical 
significance (Fig 7). Nonetheless this trend is notable because 
the surveys are deeper for later stellar types as shown in Fig 8, 
for disk radii < 100 AU. Interestingly, this trend has recently 



been found also at 10/j.m in a sample of A to M stars with ages 
between 8 Myr and 1 Gyr by Plavchan et al (2009). 

We also determine the cold disk fraction of < 10 % for 
the "old" M-dwarfs of our sample having undetermined ages 
and likely being as old as the average Galactic disk stars 
(8.8 ± 1.7 Gyr, del Peloso et al. 2005). This may indicate that 
the cold disk fraction may decrease with stellar age, similarly 
to warm disks (Rieke et al. 2005). We note that the "old" M- 
dwarfs have been observed with higher sensitivity as seen in 
Figs 5 and 6 where "young" and "old" M-dwarfs have been 
color-coded in dark and light blue. The higher sensitivity for 
the "old" M-dwarfs is because they are nearer than the "young" 
M-dwarfs. 

We now examine possible reasons why M-dwarfs have less 
detectable disks. 

Theoretical models of star formation predict that pro- 
toplanetary disks are less massive around M-dwarfs than 
around higher mass stars. Vorobyov & Basu (2008) predict 
< Md > oc < M* > 1 - 3±al , where disk mass < Md > and 
star mass < M* > are time-averaged over the star accretion 
period (0.5 to 2.5 Myr). Consequently, less primordial materi- 
als could limit planet formation around M-dwarfs. Collecting 
masses of protoplanetary disks determined by (sub)mm ob- 
servations in the nearest star forming regions, Natta, Grinin, 
& Mannings, (2000) found M d M°- 64±() I4 , but Andrews 
& Williams (2005, 2007) show that M d versus M» in Taurus- 
Auriga is so widely scattered between 0.001 and 0.2Mq that 
it precludes any meaningful correlation fit. Also, compari- 
son of disk masses of 6 members of the nearby young TW 
Hydrae Association (TWA) suggests no correlation between 
disk masses and stellar types for these reasonably coeval disks 
(Matthews, Kalas & Wyatt 2007). 

Removal of circumstellar dust by the Poynting-Robertson 
effect and radiation pressure processes are diminished around 
M-dwarfs because they are less luminous than solar-type stars, 
and so dust generated by collisions in any remnant planetesi- 
mal belt should remain there longer, giving rise to detectable 
emission. However, the opposite conclusion has been reached 
by Plavchan, Jura & Lipscy, (2005) that highlight the fact that 
dust removal around M-dwarfs could be dominated by the drag 
caused by strong winds associated with their high coronal and 
chromospheric activities. 

The formation of planetesimals and planets depends on the 
time scales between the competing processes of coagulation 
and evaporation in the early period of accumulation. Theory 
predicts that t coag increases and t evap decreases with the central 
star mass and with the strength of the FUV and EUV radia- 
tion field (eqs 47 and 48 in Adams et al. 2004, respectively). 
From their Fig 10, it can be seen that more than 10 Myr are 
needed to evaporate a protoplanetary disk around a solar-mass 
star, whereas only a few Myr are required to evaporate the same 
disk around a low-mass M-dwarf in a stellar cluster with a mod- 
erate UV flux of ~ 3000Go- This might quench planet forma- 
tion around M-dwarfs. 

Early stripping of planetesimals by passing stars is likely 
since most stars are born in clusters where stellar encounters 
as close as 160 AU are likely in the first lOOMyr (Kenyon and 
Bromley 2004b). The disruption of planetesimal disks by close 
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stellar encounters has been studied for the A6 star B Pic by 
Larwood and Kalas (2001). They found that, depending on the 
passing star's trajectory and on the relative star masses, 1 % to 
48% of the planetesimals are lost after encounters. Related is- 
sues for planetary systems were also discussed by Malmberg et 
al. (2007). 

A Lack of gaseous giant planets around M-dwarfs is pre- 
dicted by Laughlin, Bodenheimer & Adams (2004), caused by 
the longer duration required to build a core of 15M e . This lack 
of giant planets might reduce the production of second gen- 
eration dust in M-dwarf debris disks because of weaker grav- 
itationnal stirring (Wilner et al. 2002), diminishing their de- 
tectability. 

This series of arguments leads to the expectation that de- 
bris disks around M-dwarfs might be intrinsically less dusty 
and therefore more difficult to detect than those around more 
massive stars of the same age. However, most of these argu- 
ments assume that the initial protoplanetary disk mass scales 
with the central star mass, which may be plausible but is obser- 
vationally not well established. 

6. Conclusion 

To search for emission from cold debris disks, we have used the 
MAMBO-2 bolometer camera at the IRAM 30m telescope to 
map 42 nearby M-dwarfs at 1 .2 mm wavelength to a noise level 
of 0.6 to 2.8 mJy per 1 1" beam. We also reanalyzed our earlier 
MAMBO-2 and SCUBA data to form a coherent sample of 50 
M-dwarfs. Only one cold debris disk was detected, surround- 
ing the M0. 5 dwarf GJ842.2. In an attempt to discuss how de- 
tectability of cold debris disks depends on the mass of the cen- 
tral star, we have compared this result to the observed fractions 
of cold disks for more massive stars in the two submm surveys 
of Wyatt et al. (2003) and Najita & Williams (2005), who re- 
port detection rates of 22^ % for A-stars and 15^j| j % for 
FGK-stars with stellar ages between 10 and 180 Myr. For the 
19 youngest M-dwarfs (< 200 Myr) of our sample, we found a 
detection rate of 5.3^° 5 %. Hence, for this age range, there is 
a mild trend in these three detections rates, indicative of fewer 
cold debris disks detected for later stellar types - lower star 
masses - although at a low statistical significance. Nonethless 
this trend is notable because the sensitivities of these surveys 
are deeper for later stellar types. We also determine the cold 
disk fraction of < 10 % for the "old" M-dwarfs (likely a few 
Gyr) of our sample, indicative that the cold disk fraction may 
decrease with stellar age, as is also seen for warm disks. Future 
observations of a larger and better controlled sample of stars 
of all stellar types with Herschel in the far-IR and deeper ob- 
servations in the (sub)mm will be able to better clarify these 
issues. 

7. Appendix : computation of the dust temperature 
and mass for the possible large debris disk 
around GJ526 

The source that heats dust in debris disks is usually the stel- 
lar radiation field, but in a large debris disk, as possibly found 
around GJ526, the interstellar radiation field dominates at some 
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Fig Al : Variation of the temperature of a grain 100//m in diameter 
exposed to the radiation of the stellar field and of the isotropic inter- 
stellar radiation field. This latter field dominates at some disk radius R 
that depends on the stellar luminosity. The saturation temperation of 
4.9K due to the interstellar field is consistent with the computation of 
Kriigel (2003, p. 249). For M6 only, the dashed line indicates the grain 
temperature when the interstellar radiation field is not included. (This 
figure is available in color in electronic form). 

radius from the star. We compute this transition radius for sev- 
eral stellar spectral types by solving numerically the integral 
equation for a grain at thermal equilibrium absorbing both the 
stellar and interstellar incident fields : 

X+oo 
Qab.M, a) . (na 2 J t (A, r)+4na 2 J ISRF (A)) dA = 

r>+oo 

| 4na 2 Q abs (A, o)ttB(A, T g )dA (Al) 
Jo 

where Q a bs(A, a) = 1 for A < Aq and Q a bs(A, a) = Ao/A with 
Ao = 2na for A > Ao, which approximate the absorption ef- 
ficiency computed for carbon by Laor & Draine (1993). The 
parameter a is the radius of a spherical grain. The expression 
B(A, T g ) is the Planck function at grain temperature T g that de- 
pends on a. The intensity of the stellar radiation field J,(A, r) is 
nB(T t , A) x (AnRt/Ajiry at disk radius r with the star character- 
ized by its effective temperature T» and its radius R„ . Intensity 
of the Interstellar Radiation Field Jisrf(A) m the solar neigh- 
bourhood is : 

i=6 

JisrfW = Yj CinB{A, 27) 

i=0 

with components C = l,7o = 2.1K,C\ = 4x 10~ 17 ,:Ti = 

27000/:, c 2 = io~ 14 ,r 2 = isook, c 3 = io~ 13 ,r 3 = 
40oo^,c 4 = 4x io- 13 ,r 4 = 3000^,c 5 = 1.5 x io- 1() ,r 5 = 

400K, C 6 = 0.5 x 10~ 5 ,r 6 = 46K, at galactocentric distance 
lOkpc (Mathis, Mezger & Panagia 1983). In Fig Al, we solved 
numerically eq (Al) to determine dust temperature T g (r) as 
a function of disk radius in modelling grains with the single 
size 2a = lOO/vm, typical for submm observations, and for 
several stellar spectral types. For GJ526 (Ml. 5), this figure 
shows that the interstellar radiation field becomes dominant at 
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Fig A2 : Variation of temperature as a function of grain radius at 
500 AU from the Ml. 5 dwarf GJ526. Similar results are reported by 
Kriigel (2003, p. 249). The discontinuity in the silicate curve arises 
from the piecewise-defined funtion Q a bs(A,d) adopted from Laor & 
Draine (1993). (This figure is available in color in electronic form). 



r > 1000 AU, but this transition radius is only 200 AU for 
spectral type M6. 

In Fig A2, we use a more realistic model for the dust in 
adopting the standard grain size distribution dN = Noa~ 3 5 da 
(Dohnanyi 1969) to compute temperature T g (a) as a function 
of grain size at r — 500 AU with eq (Al) but modified to 
include the grain size distribution. Temperature increases sig- 
nificantly for a < 10 /im which is relevant for M-dwarfs be- 
cause the grain blow-out size is small. For GJ526 (0.03 1L ), 
it is as small as 0.04 /jm for carboneous grains, ignoring the 
effect of stellar wind drag for this star which has a low coro- 
nal/chromospheric activity (LogLj=26.87 [10~ 7 W], Schmitt & 
Liefke, 2004). Finally, we use this temperature function T g (a) 
at r — 500 AU to determine the dust mass around GJ526 by 
matching the total flux density of the 5 clumps (S 1.2 mm = 
21.2 + 2 mJy) to the flux density of our model. If d is the dis- 
tance to the star, the predicted flux density is : 



No 
And 2 



r 



Ana 2 Qabs(v, . nB(v, T g (a)) . a 3-5 da 



(AT) 



the limit «„„•„ is set by the blow-out size (a mi „ = 0.04/im 
for GJ526). The limit a max is related to the total dust mass 
Md - yja max for the size distribution above and spheri- 

cal grains of density p. Md is the dust mass probed by the mea- 
sured flux density at the observing A. Practically, the grain size 
limit a max is the value that makes convergent the computation 
of the power emitted by Md over a band 2b centered on the 
observed wavelength (1.2 mm). We adopted the convergence 
criterium of 5% to match the relative accuracy of the measured 
flux density. In other words, a max is increased until the integral 
below converges to within 5% : 



Fig A3 SED of carboneous dust around GJ526 based on our model fit- 
ted to the 1.2 mm flux density. Our model includes the collisional dust 
size distribution dN = a da and the resulting non-uniform temper- 
ature for grains of various sizes computed in Fig A2. Superimposed is 
the corresponding grey-body model computed with the single grain 
temperature 4.9 K and standard mass opacity 1.7 cm 2 g~' at A = 
850 [im scaled with the efficiency 210 fim/A (see text). (This figure 
is available in color in electronic form). 



only if the flux density were measured more accurately. We 
found that a max is ~ 24mm and ~ 36.4mm for carboneous and 
silicate grains, respectively, and, independently of the choice 
b = 5%, 10%, 20% x 1.2 mm. Similar calculation was made by 
Wyatt & Dent (2002) for the A3V dwarf Fomalhaut by integrat- 
ing not only over grain size but also over a range of radius r to 
account for spatial distribution of the dust in a ring. They found 
that 95% of the flux density comes from grains and pebbles less 
than 100mm in radius. This is comparable to our determination 
of a max and the difference must come from different function 
for T g (a) and their hypothesis of dust spatially distributed. 

Figs A3 and A4 show the SED of the dust around GJ526 
based on our model for both carboneous and silicate grains 
adjusted to the flux density measured at 1.2 mm. The dust 
mass determined by our model is between ~ 6 and ~ 10 lu- 
nar masses. We have added in these figures, the SED based on 
the standard grey-body model of optically thin dust emission 
(Hildebrand 1983 and Zuckerman 2001) : 

S v = M d x B(v,T„) x Kabs/d 2 (A4) 



where the mass opacity K a t s is conventionally taken to be 
1.7 cm 2 g~' at 850 yum and is scaled by 210 fim/A for A > 
210 fim. This grey dust mass is about 2 to 4 times larger than 
the one determined by our model. This level of accuracy is typ- 
ical of dust mass estimates for debris disks, presently. 

The SED of our model is significantly more extended to the 
Far-IR and to the radio than the grey-body model when both 
models match the flux density at A = 1.2 mm. Noticeably, we 
point out that this Far-IR extension is much less pronounced 
~ 3 5 da dA (A3)if our model is applied to an A-star. Future observations with 
Herschel will test these SED differences and thus will probe, in 
This power would need to be computed more accurately by in- essence, the collisional dust size distribution adopted and the 
eluding the small contribution to the emission of larger pebbles non-uniform grain temperature. 



X+1.2 mm+b na max 
I Ana 2 Qabsi^j a) ■ 7rB(A, T g (a)) . a 
1.2 mm-b J a min 
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Fig A4 SED of silicate dust around GJ526 based on our model and 
fitted to the 1.2 mm flux density. Same comments as in Fig A3. (This 
figure is available in color in electronic form). 
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